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Abstract 

We reconsider string and domain wall central charges in M = 2 super- 
symmetric gauge theories in four dimensions in presence of the Omega 
background in the Nekrasov-Shatashvili (NS) limit. Existence of these 
charges entails presence of the corresponding topological defects in the 
theory - vortices and domain walls. In spirit of the 4d/2d duality we 
discuss the worldsheet low energy effective theory living on the BPS vor- 
tex in J\f = 2 Supersymmetric Quantum Chromodynamics (SQCD). We 
discuss some aspects of the brane realization of the dualities between 
various quantum integrable models. A chain of such dualities enables 
us to check the AGT correspondence in the NS limit. 
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1 Introduction 



Recent work by Nekrasov and Shatashvili [T] has initiated the program of quantization of 
integrable systems by deforming four dimensional supersymmetric (SUSY) theories these in- 
tegrable systems are associated with. The relationship between classical integrable systems 
and N = 2 supersymmetric gauge theories has been extensively studied for a long time (2j|3). 
Remarkably the low energy dynamics of M = 2 gauge theories is captured by finite dimen- 
sional integrable systems such that the phase space of the latter is related to the instanton 
moduli space of the former. To be precise, there are two different classical integrable systems 
involved here. The first one is a holomorphic integrable system of the Hitchin or the spin 
chain type giving a Seiberg-Witten curve whose Jacobian is mapped onto complex Liouville 
tori and action variables are identified with order parameters in the gauge theory. The sec- 
ond integrable system is of the Whitham type emerges from the renormalization group (RG) 
flows, and the very RG equation plays the role of the Hamilton- Jacobi equation written in 
proper variables. 

The quantum integrable systems we are considering in this paper can be extracted from 
four dimensional theories in Omega background [4| with ei = e, £2 = (the so called 
Nekrasov-Shatashvili (NS) limit [5]). Given a prepotential of the 4d theory J-"(a, ex, £2) as a 
function of the Coulomb branch moduli parameters {a} and the Omega deformation param- 
eters ex t 2, we can consider in the NS limit an effective 2d theory with the following effective 
exact twisted superpotential 

W(a,e) = hm ^ 6; 62) , (1.1) 

where e\ is replaced by e |1|. The F-terms of the effective Lagrangian effectively become two 
dimensional in the NS limit and are described by W(a,e). For small e formula (1.1) can be 
even further simplified 

W(a,e) = ^ + ..., (1.2) 

where the ellipses denote terms which are regular in e. The above twisted superpotential 
includes both perturbative and instantonic contributions. 



Minimization of superpotential (1.1) yields the supersymmetric vacua which, according 
to the same authors [6j[5] are intimately connected with quantum integrable systems. Indeed, 
according to Nekrasov and Shatashvili, supersymmetric vacua of an appropriate two (also 
three and four) dimensional M = 2 gauge theory are in one-to-one correspondence with the 
Bethe roots of a certain integrable_system. It is useful to make the Legendre transform and 
consider the dual superpotential Wd(cid,^) depending on the dual variable ap- Thus the 
equation 

ex P ^ =1 > ( L3 ) 



da 



D 



can be viewed as a Bethe ansatz equation for some integrable system. This is consistent 
with the interpretation of the prepotential as an action in the Whitham system. Since ap is 



the coordinate variable in the Whitham dynamics, one can recognize the exponent in (1.3) 
as the canonical conjugate momentum. 
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There also exists a different well known duality between four dimensional gauge theories 



and two dimensional gauge theories (linear sigma models) [7}|9] (see 10 for review, we shall 
refer to it as 4d/2d duality). The four dimensional theory here sits at the root of its baryonic 
Higgs branch; therefore electric and flavor charges can be combined in a single set of quantum 
numbers. Together with magnetic charges they form two sets of quantum numbers which 
parameterize masses of four dimensional dyons. On the two dimensional side one has kinks 
interpolating between different supersymmetric vacua, their masses also depend on two sets 
of charges - Noether and topological ones. The statement of the 4d/2d duality in its original 
formulation [7] is that the two sets perfectly match with each other. Thus the BPS spectrum 
of an entire four dimensional theory can be studied using a relatively simple two dimensional 
(gauged linear) sigma model with four supercharges. It was later shown by Shifman and 
Yung [To] that the 4d/2d correspondence is not accidental, the underlying two-dimensional 
theory in fact should be treated as a low energy effective theory on the worldsheet of a BPS 
vortex. 

The first step in the direction of merging the 4d/2d correspondence and the 
gauge/quantum integrabilty duality together was put forward by Dorey, Hollowood and 
Lee [11 1 . The authors considered all three ingredients at once - four dimensional J\f = 2 
gauge theory (SQCD with Nf = 2N), a two dimensional Af = (2,2) gauge linear sigma 
model (GLSM) (a certain U(K) gauge theory), and an integrable system. Critical ingredi- 
ent of the duality is that the 4d gauge theory sits at a baryonic root of the Higgs branch, 
which in Omega background undergoes a deformation, and scalar field VEV gets shifted by 
an amount proportional to e. Provided the baryonic Higgs root condition is satisfied, the 4d 
theory is shown to be dual to a given GLSM whose twisted superpotential plays a role of the 
Yang- Yang function for a SL(2, R) Heisenberg magnet. The four dimensional superpotential 
is shown to be equal to the effective 2d twisted superpotential of the form (1.1) on shell. 

The initial motivation for this research was the identification of the different dualities 
known for quantum integrable systems in the framework of SUSY gauge theories and their 
brane realizations. However, first it is necessary to explain the geometrical meaning of 
degrees of freedom in relevant integrable systems which is a subtle issue. It was clear for 
a while that these degrees of freedom are described in terms of brane embeddings into the 



internal space. We refer the reader to 12 , where geometrical aspects of dualities between 
integrable systems are reviewed. 

A proper brane content involves surface operators, or equivalently nonabelian strings 
with large tension. To begin with we revisit the classification of the BPS solitons of different 
codimensions in the e deformed theory. There is some surprise. It turns out that some 
BPS states in the deformed theory to the best of our knowledge were overlooked. We study 
carefully the central charges of the Omega-deformed SUSY algebra and argue that there 
are new stringy and domain wall type central charges. The key point is that tensions of 
such strings and domain walls are proportional to the graviphoton field and these defects are 
absent in the undeformed theory. We shall investigate the string-like object both for the pure 
supersymmetric Yang-Mills (SYM) theory and for the SQCD. It will be argued that there 
are some singularities akin to those of cosmic strings. The corresponding BPS equations for 
such strings will be derived and finiteness of their tension will be discussed. Similarly the 
domain walls solitons will be found, which are more expected to appear, since due to the 
Omega deformation the theory has the discrete set of vacuum states, hence domain walls 
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naturally emerge. The monopoles in Omega background which have been already discussed 
in the literature 13, 14 also require some attention for their proper interpretation. 

Turning to the dualities in quantum integrable systems we shall focus on two subjects. 
First we shall consider the quantum bispectral duality relating two different integrable sys- 
tems. Classically eigenvalues of the Lax operator in one system get interchanged with coor- 
dinates in the second system. Quantum mechanically it means that the single wave function 
serves for two systems simultaneously when considered as the function of the spectral or 
coordinate variables. Although this question has not yet been elaborated to the full extent 
in the literature, by employing the quantum version of the duality we were able to explain 
the details of the 4d/2d duality in Omega background [11] . Geometrically bispectrality cor- 
responds to a rotation of the brane configuration which represents the 4d gauge theory in 
question. We consider the relation between Bethe ansatz equations (BAE) for dual inte- 
grable systems and briefly discuss their degenerate solutions corresponding to the analogues 
of Argyres- Douglas points. Using the duality between the families of spin chain models and 
the Calogero-Ruijsenaars systems we shall identify bispectral pairs in both families. As a 



byproduct we will show that the Alday-Gaiotto-Tachikawa (AGT) duality 15 in the NS 
limit can be proved using the chain of dualities involving the bispectral pair. 

The paper is organized as follows. In the next section we review how the M = 2 su- 
persymmetry algebra is affected by the Omega background. We then compute the central 
charge for a BPS string in pure M = 2 Super Yang Mills theory and calculate its tension. In 
Sec. [3] we consider domain walls and monopoles in Omega deformed SYM. Then in Sec. [4] we 
investigate BPS strings in well-studied example of the 4d/2d duality - the supersymmetric 
QCD. Section [5] is devoted to brane constructions and integrable systems associated with 
M = 2 gauge theories and dualities between them. In Sec. [6] we show that in the NS limit, 
the celebrated AGT correspondence 15 can be reduced to the so-called bispectral duality 
between two integrable systems. Finally in Sec. [7] we conclude and speculate on further 
research topics. 



2 Flux Tubes in Pure Super Yang-Mills Theory 

In the standard lore of topological defects in supersymmetric theories, the BPS strings only 
exist when a gauge group is at least semi-simple, e.g. U(N). A simple reason for this is 
based on existence of a nontrivial fundamental group of the resulting moduli space due to 
presence of a U(l) factor. The latter causes a nonzero Fayet-Iliopoulos (FI) term which 
supports string solutions, and we shall refer them as FI strings. In this section, we shall 
instead consider to a new kind of string-like objects which have not been discussed in the 
literature before, we shall refer to them as e-strings. As we shall later see their tension is 
proportional to e 2 and construct the relevant classical field configurations. For simplicity we 
shall only focus on the gauge group SU (2) in this section. 

Action. Let us start with the M = 2 Super Yang-Mills theory in four dimensions in Omega 
background. To set the notations, the Lagrangian of the undeformed theory reads 



Tr J d 4 6$e v <$>e- v + Tr J d 2 9 (W a ) 



(2.1) 
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where $ = (0, if), F) is adjoint chiral superfield, V = (a, A, D) is adjoint vector superfield, 

Aizi I e 
9 2 ^ 2tt 



W a is its field strength, and r = ^? + is coupling constant. 



Omega deformation. The Omega deformation of a four dimensional theory like (2.1 ) can 
be constructed from a six dimensional theory by compactifying the theory on a two-torus 
with twisted boundary conditions 16 17 . Torus action on IR 4 is given by two matrices Q™ n 



where m, n = 1, 2, 3, 4 and a = 5, 6 which act by rotations in 12 and 34 planes respectively. In 
the NS limit matrix Q$ vanishes, therefore we shall denote Q = f2 5 . Metric on the deformed 
torus reads 

G AB dx A dx B = Adzdz + {dx m + n m dz + tTdz) 2 , (2.2) 

where z = x 5 + ix 6 , z = x 5 — ix e and the vector field Vt m = Vt™x n . In the notations of [l3] 
Qm _ ^_i ex 2^ i ex ^^ q ? q)_ j n other words vector field Q = ied v is a rotation generator around 

X3-axis. Here we denote p = yx\ + x\. The components of the metric in the limit A — > 
read 

Upon the dimension reduction, the fifth and sixth components of the gauge field form an 
adjoint scalar, which undergoes the following deformation due to the Omega background 

_> _ iQm Vm + Ltfnns^ , ( 2 . 4 ) 



where 0, m , Vt mn were introduced after formula (2.2) and S mn is the spin operator for adjoint 
representation of the gauge group. The latter does not affect the bosonic part of the theory, 
however it does modify the fermionic part. This issue will be important when we consider 



the supersymmetry algebra of the theory momentarily. Transformation (2.4) itself is not a 



well-defined change of coordinates, but <f> enters the Lagrangian in a special way, this shift 



brings us to a well defined Lagrangian of a modified theory 18 . Another deformation of the 
theory consists of shifting of the coupling constant, thereby we promote it to a superfield. 
In the M = 2 superfield language^] the shift reads as follows 

r^r-e m e n (a mn )\ (2.5) 

where 9 m = (& m ) aI 6ai is the twisted Grassmann variable for the diagonal su(2)^ + 7£ gener- 
ators. In components the Lagrangian of the M = 2 SYM after the deformation takes the 
following form 



4- > .... .~ . 

where / = 1, 2 denotes the R-symmetry index, and spinor indices are suppressed. 



) s \f a (Q m V m - ln mn a mn )X a f - y\ a f {n m V m - \9r n a mn )V a , (2.6) 



See Shadchin's PhD thesis 17 for details 
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SUSY transformations. Recall that M = 2 super symmetry algebra in four dimensions 
has the following form 



{Qai Qja} — 2Paa$J + 2Z a aSj , 
{Qai Qp} = e a/3£ IJ Z mon + (^d.w.)a^ ' (2-7) 

There are three types on central charges: string, monopole and domain wall types. We shall 
focus on the former in this section leaving monopoles and domain walls to Sec. |3j 

The full global symmetry of the theory is SU(2) L x SU(2) R x SU{2) n (left, right and the 
R-symmetry). It is broken by the Omega background in the NS limit to SU{2) L x SU(2) R+n 
by paring the R-symmetry with the right handed SU(2). The supercharges undergo the 



Donaldson- Witten twist 19 



Q — $lQa ' Qm — {&m) Qla > Qrnn — {&mn)lQa • (2-8) 

These transformations can be inverted as follows 

Qa = \{ am )aQm ^ Q 'aj = \^ajQ + \{o mrl ) ajQ mn ■ (2.9) 

It turns out that a generic Omega background breaks all supersymmetries of the theory 



(2.6) except the BRST charge Q. Moreover, it can be shown that the Lagrangian (2.6) is 
a Q-exact expression which makes it possible to compute the partition function of the 
theory by localization methods. 

It is more or less clear that the obstacle to supersymmetry is due to the spin operator 
terms ^Q mn a mn in the fermionic sector. The theory thus has to be further deformed to gain 
more supersymmetry. To understand what we need to do, let us look at the supersymmetry 



transformations of (2.6 ) with the problematic spin operators omitted, and see what additional 



terms do we need to introduce. One has 13 the following under variations 



+ C //3 (O^(V m - F mn Q n ) . (2.10) 
However, since the action is not supersymmetric, the above transformations do not leave the 



Lagrangian invariant. As suggested in 16 and later extended in 13 , one has to turn on 
R-symmetry Wilson lines properly to restore partial supersymmetries. Thus, in the NS limit 
one has to add 

-AjX'Xj-Aj^Xj, (2.11) 

where 

A J j = -^ mn (a mn yj , (2.12) 

to the Lagrangian. One can also treat the above terms as emerging from a superpotential 
which we add to the theory, and we shall speculate more on this in Sec. [3] In the NS limit, 
when 62 = 0, the supersymmetry of the theory is enhanced to J\f = (2, 2) [13] and is generated 
by the following supercharges 

Qi,Q2,Qi3,Qu- (2.13) 
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Using the inverted transformation (B.3 ) we conclude that in the original formulation of SUSY 
algebra (B.l ) the following generators are included into Af = (2, 2) subalgebra 



Q12, Q21, Q i2 , Qn- (2.14) 

In the remaining part of the section we shall investigate 1/2 BPS object - a string which is 
annihilated by the above charges. 

String central charge and string tension. The supercurrent for the Omega deformed 
SYM theory was computed in [13]. Its Euclidean time component has the following form 
(assuming static configuration, B 3 7^ 0, others components of F mn vanish) 

+ ^(^ 4n )«(-Vn0 + F np QP)\ ip (2.15) 

Let us find the string central charge current. Performing standard variation of the above 
supercurrent we getn 



(2.16) 



where P4 is Hamiltonian of the system. Note that there is an additional contribution to the 
above string charge current which is bilinear in fermions of the form d m (Q rn \\). For classical 
analysis, where all fermionic fields can be put to zero, this contribution can be omitted. We 
see that there is a correction which represents the string central charge, specifically the 
correction takes the following form 



where p 2 = x\ + x\ is the transversal coordinate to the string. If e is real then 



c 3 = M>(9W° s 3) 



The central charge is given by 



2rr 



J string 



d 3 x( 3 = j 2 jdzjdppj dipd v {iRe(e<j> a )B*) 
dz I dppB%Me(e(, 





27T 



2 Technically there is another contribution from the R-current 
charge.. 



20 



which contributes to the ( 



(2.17) 



(2.18) 



1 1 



(2.19) 

central 
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We can immediately see that multi-valuedness of as a function of the azimuthal angle 
is required in order to make the central charge nonzero. The tension of the string solution 
under consideration (let's call them e-strings) is therefore given by 



Assuming that 



where a is an constant, we arrive to 



dppme (eB%<j) a (e- 2wia - 1)) 



(2.20) 



(2-21) 



(2.22) 



The above expression for the tension of e-string only makes sense if it is finite. In order to 
establish that one has to solve BPS equations in order to find the profile functions for and 
B 3 as function of the radial coordinate p. 



BPS equations. Let us now find the BPS equations which describe such a string. Once 
supersymmetry algebra is understood (2.14), we can focus on the bosonic part of the action 

i|2 



4g 2 mn 



V m 6 - F, 



Q"|2 i J_ 



- iV m (ft r 



(2.23) 



Note that in the NS limit Q m Q n F£ in identically vanishes. Let us now do the Bogomolny 
completion, as the supersymmetry algebra suggests 



C 



i 

2<? 2 



I 2 + o^|Vi0 a + iV 2 <p a - (fi 2 - itii)B. 



a 1 2 
3 I 



2g 2 



(2.24) 

The above inequality is saturated provided that the following BPS equations are satisfied 

Bl + (f)T a (f) - zV m (fT> a - tt m (f) a ) = , 

Vi0 a + iV 2 <f) a - (0 2 - zfii)B£ = . (2.25) 
One can also check that the above BPS equations are consistent with the M = (2,2) super- 



symmetry algebra. Indeed, by looking at the gluino variation in (2.10) we need to set to zero 



all expressions which enter the right hand side together with £n, (22 and their complex con- 
jugates. Contributions proportional to £12 , £21 and conjugated terms vanish automatically 



due to the BPS condition. By doing so one arrives at equations (2.25) 



Sometimes it is more convenient to switch to the complex coordinates 

x + iy , w = x — iy , 



w 



than the BPS equations (2.25) take the following form 

d^Al - d w A% + e a V0 c - j(wV, - wV w )(e4> a - e<j> a ) 
d _ r _ e^AlP + ewid^Al - d w A%) 








(2.26) 



(2.27) 



The above equations can be used in study of the effective two dimensional theory living on 
the e-string. 
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Solution of BPS equations and vortex tension. Let us proceed with the solution of 



(2.25). We will look for a background solution when all fields are aligned along the Cartan 
subalgebra of the gauge algebra. Thus color superscript will always be a = 3, in the rest of 
the section we shall omit it. One has the following 

B 3 - id m {Q m - tt m <f)) = 0, 
<9i0 + id 2 (j) - (fi 2 - ifii)S 3 = . (2.28) 

Decomposing <fi = (pi + z'0 2) e = ei + ze 2 into real and imaginary parts we obtain 

B 3 = 29 v (ei0i + e 2 2 ) , (2.29) 

and two first order equations on <f>i and 2 . After some simple manipulations one gets 

A0i + 4e 2 ^(ei0i + e 2 2 ) = , 

A0 2 - 4ei^(ei0i + e 2 2 ) = , (2.30) 

which, after adding these equations with proper coefficients, implies that e\<f>x + e 2 2 is 
a harmonic function. However, at this point boundary conditions of the solution remain 



unclear as we need to have 0(27r) ^ 0(0) in order to gain finite tension (2.22). In order 
to make the problem mathematically precise we can make the following trick. The phase 
difference of e 2ma will be identified with the deficit angle of a cone which is obtained by 
gluing ip = ray with <p = 2tt one. 

Solution of Laplace equation on a cone is formally given by a series of positive and 
negative powers of p with angle dependent coefficients. The latter are normally expressed in 
terms of ellipsoidal harmonics. Since we are interested in normalizable solutions we only leave 
negative powers of the radial coordinate in the series. The solution will however be divergent 
at the origin. The dependence on phase a is then hidden in the harmonic coefficients. We 
shall refrain from giving more details here since we are not investigating any dynamics on 
e-strings in this paper. Rather we provide the evidence of existence of e-strings and finiteness 
of their tension. 



Let us now look at the vortex tension (2.22). Using (2.29) we conclude that 



2-k 

1 f d 2 xTiB 2 = ^ J dip J dppd^e(^ a )B^. (2.31) 



2g 



As we argued above the integral over the radial coordinate diverges, however it does not 
make the tension infinite. To see this let us regularize the radial integral on its lower limit 
by putting a cutoff at some small value of p = po- After integrating the full angular derivative 
we see that the contribution to the integral at the lower limit cancel each other as po — >* 0- 
Thus we are left with the contribution from large p. As we are not specifying the full solution 



of the BPS equations we shall not evaluate integral (2.31) here. However, from dimensional 
ground we anticipate 

T = A(a)\e\ 2 , (2.32) 

where A is a constant. Therefore we find that the tension of e-strings is quadratic in e, and 
it is only nonzero for fractional winding numbers. 
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To conclude this section let us make a few general comments concerning e-strings. First, 
at large values of the graviphoton field the string tension is large and we can safely consider it 
as semiclassical object and the string could serve as the new type of the surface operator. In 
this case one can use the standard technique to get the worldvolume theory. We shall discuss 
the worldvolume theory elsewhere. The second point to be mentioned is some analogue with 



the string in the noncommutative gauge theory found by Gross and Nekrasov 21 . They have 



discovered that the Dirac string attached to the monopole in the noncommutative theory 
becomes observable and its tension is proportional to the noncommutativity parameter. Since 
using the chain of dualities [22] the e parameter can be traded to the noncommutativity in 
the internal space one could look for more close relation between two types of strings. Finally, 
one could ask for the brane realization of the e-string. Certainly it can not be identified as 
D2 brane similar to the FI string since we can not reproduce the tension with such brane 
realization. Hence the most natural candidate is the properly embedded D4 brane. We hope 
to discuss the details of the brane realization of e-strings elsewhere. 



3 Monopoles and Domain Walls 



We have discussed color flux tubes in Omega background in Sec. |2j Here we shall address two 
other types of topological defects we often encounter in supersymmetric theories - monopoles 
and domain walls. 



Central charge. 



Recall from (B.l ) that domain walls and monopoles saturate holomorphic 

Symmetric combination of these charges give 



central charges in supersymmetry algebra, 
domain wall piece, whereas an antisymmetric one contributes to monopoles. Analogously to 



a string charge density we have evaluated in (2.16), we can proceed with the monopole and 
domain wall. The former calculation has been performed in 13 , and the latter gives 



& bp J Jo 



(3.1) 



We have included here only the bosonic contribution to the supercurrent which is relevant for 
classical analysis. The full expression will also contain bilinear term in fermions d m (Q rn XX) + 
H.c. which will be manifest for quantum calculations [23]. 



Structurally (|3.l|) is very reminiscent of the string current (2.16). As in the string charge 

") , (3-2) 



Case, cL field dependent FI term appears 
leading to the following expression of the domain wall tension 

rp 1 C a ( A a 



(3.3) 



which can be viewed as a non-Abelian generalization of the standard calculation in theories 
with superpotentials. Let us now have a closer look on the BPS monopoles and domain walls 
appearing in M = 2 SYM theory, again for simplicity we shall consider 577(2) gauge group 
here. 
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BPS equations and monopole solution. Ito et al [13, 14 have investigated BPS 
monopole solution of N = 2 SYM with gauge group SU (2) in Omega background in the NS 
limit. It reads 

Bf- V i( f> a + ie ijk WB ak = 0, (3.4) 

or in components 

#3 -V 3 a - ex m B^ = 0, 

B a m -V m <p a + ex m Bl = Q } (3.5) 



where m = 1, 2. The solution of these equations is given in 13 . The authors' conclusion is 
that the monopole's mass is not changed, the magnetic field strength has the same form as 
the one in the undeformed case for e = 0. However, there is a correction of the scalar field 
profile. In the singular gauge (when only <p 3 ^ 0) the solution reads 



3 (p, z) = K \v 2 H{p, zf + e 2 + 2 ^=^$ + G{p, zf , (3.6) 



z 2 



where function 

2p 2 vzeF(p,z)H(p,z) p 2 z 2 e 2 F(p, z f 2p 2 e 2 F(p,z) 



G(p,z) 



p2 _|_ z l pi, _|_ z l pi _|_ z l 

_ 2p 2 vzeF(p,z)H(p,z) p 4 e 2 F(p,zf 

(p 2 + z 2 f 2 (p 2 + z 2 f [ } 

vanishes at z — > ±oo. Functions F and H are taken from the 't-Hooft-Polyakov monopole 
solution |24| 



H(p,z) = coth^y^T^+ , F(p,z) = l (3.8) 

vyp 2 + z 2 sinh t> a/p 2 + -2 2 

Note that scalar field does not go to its vacuum value v any longer as it does for e = 0, 
but rather interpolates between <f) +OQ — v + e to 0_oo = t> — e at plus and minus z-infinity 
respectively. This suggests us that maybe 1/2 BPS monopole is not a proper interpretation 
of the above solution and more structures can be involved. 

Before we go further, let us mention an useful symmetry of equations (3.5). In the above 
analysis axial symmetry was assumed such that both cf) and B fields depended only on p and 
z. We can also introduce azimuthal angle <p in the game by giving the scalar field a phase 

(j)^ (j)e ia<p . (3.9) 



In order to preserve (3.5) the magnetic field strength also acquires a phase and its azimuthal 
component B® gets generated due to V m term in the second equation. Provided that such 
a configuration is chosen, the FI field reads 

C = f 2 a^ a + 4 a )- (3.io) 
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Then the corresponding domain wall's tension (3.3) becomes 



T 



r eTr0 2 



(3.1i; 



for imaginary and real e. Since (3.6) still solves equations (3.5) we can compute the tension 



for the solution in hand. One gets 

T=^e{{v + ef-{v-ef) = Ave* 



(3.12) 



Let us now see how to construct monopoles and domain walls in SU(2) SYM theory in the 
NS Omega background. 

Monopole on a domain wall. We shall perform a slightly different Bogomolny comple- 



tion of the action (2.23) than the authors of 13[ 14 
+ ^\B a 1 +iB a 2 + {V 1 + iV 2 )cj ) 



ak\ 2 



(n 2 -itti)B%\ 



+ ±d m (B*(n m $ a - n m <f> a )) + fM<P a ^ m {^ m r - n m r)) - ^(W) . (3.13) 

Three terms in the third line above correspond to strings, domain walls and monopoles 
respectively. Existence of the first two types of solitons solely relies on the non-trivial field 



dependent FI term (3.2) 



(3.14) 



BPS equations follow from (3.13) immediately 



(3.15) 



Bl + V 3 a + ie 3jk tt j B ak - iV m (tt m <p a - fl m (f) a ) = , 
B{ + iBl + (Vi + iV 2 )0 a - (Q 2 - = . 

We can now observe that if e and <p a are real valued, as they were chosen to be in [l3|[l4] , 
then the latter equation above splits into two, one for the real part, one for the imaginary 
part of its l.h.s. Also the field dependent FI term vanishes. We immediately identify them 



as the last two equations of (3.5). However, if a different ansatz is chosen, when either e or 



i a or both have imaginary part, the FI parameter (3.2) kicks in, and equations (3.15) no 



longer decouple. Their solution is probably more complicated than reported in 14 , and will 
be reported elsewhere. 

Boojums. At this point let us also mention that the setup we have just described also 



admits strings provided that the parameter a in (3.9) is non integer, otherwise the string 



central charge vanishes, indeed it follows from (2.22). Thus the solution above describes 



a BPS monopole on a domain wall, and, if a is not integer, a more complicated boojum 
10 which is a junction of a string, domain wall and a monopole Fig. [I] it is a 1/4 



construction 



BPS configuration. Remarkably, all three structures coexist together and emerge together 
from Omega deformation. It appears to be impossible, as far as our analysis suggests, to 
find, say, only domain wall without a monopole, or vice versa - they always come in pairs. 
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Figure 1: Left: Boojum as a monopole-string-domain wall junction. The string is infinite and is 
stretched along the z-axis. Right: e-string ending at monopoles located on two parallel domain 
walls in xy-plane can be viewed as the superposition of two boojums. The string does not continue 
through the domain walls to the outer area, since the scalar field, main building block of the e-string, 
vanishes outside of the domain walls. 



Strings, however, as we discussed in Sec. |2j can exist on their on provided that a is a non 
integer. 

Interestingly, the object we have just described - a boojum, can be placed on ends of e- 
strings, the same way as monopoles in the Higgs phase oiN = 2 theory can have non-Abelian 
flux tubes emerging from them [To] . Such a string is depicted in Fig. [T] 



Relationship to coupled 4d/2d systems. In 25,26 coupled 4d/2d systems were stud- 
ied. An example of such a system is a four dimensional gauge theory with a surface operator 
insertion. The 4d theory is considered to be in the Coulomb branch, a 2d theory lives 
on the surface defect and both systems are coupled. Remarkably, both 2d (kinks) and 4d 
(monopoles, dyons) BPS states can be found in such systems and the authors of [26] managed 
to derive the full 2d/4d wall crossing formula. Bound states of monopoles on surface defects 
are present in the theory, since the 4d theory is at the Coulomb branch, its magnetic field 
has a spherically symmetric pattern, unlike a Higgs monopole whose field lines are trapped 
to a vortex. These two pictures - Higgsed monopole an a vortex and a Polyakov-'t-Hooft 
monopole on a surface defect Fig. [2] may represent two different limiting configurations of a 
more generic setup, which involves more sophisticated 2d/4d dynamics. Keeping the calcu- 
lations we have done in this section, we may hope that 4d theories in Omega background 
may be reasonable candidates for such a theory. It would be interesting to investigate the 



solution of BPS equations (3.15) more closely and study different values of the deformation 
parameter e. At large e the surface operator limit emerges and Gaiotto et al story |26| may 
also arise. 
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Figure 2: 't-Hooft-Polyakov monopole on a surface defect. 

4 Af = 2 SQCD in Omega Background 

The 4d/2d duality was initially formulated for J\f = 2 supersymmetric QCD [7], with gauge 
group U(N) with iV < iVy < 2N flavors. Physical explanation of the duality [9j[8j (see 
also [27]) relies on existence of BPS flux tubes (FI strings) - solitonic solutions of the theory. 
Below we demonstrate how these solutions can be constructed and what are the resulting 
BPS equations for Omega-deformed theory, which as we know, preserves only M = (2, 2) 
part of supersymmetry (4 out of 8 supercharges). Similarly to the pure SYM discussed in 
the previous section, the BPS vortex in question will invariant under the unbroken part of 
the SUSY algebra, therefore it will still remain to be a BPS configuration. 



4.1 Chen-Dorey-Hollowood-Lee duality 

Here we shall review main aspects of the duality unveiled in [11,28 . The authors have 
proposed and proved that in NS limit, four dimensional U(N) SQCD with N fundamental 
hypermultiplets of masses mi, ...m^ together with N antifundamental hypermultiplets of 
masses rhi, . . . fh^ and coupling constant r is dual to the two dimensional U(K) GLSM with 
N chiral fundamentals of twisted masses Mi, . . . , M n together with N chiral antifundamentals 
of twisted masses Mi, . . . , M n and coupling constant f. The above statement holds provided 
that the four dimensional theory is considered in the Higgs branch defined by the condition 

<p a = m a -n a e, (4.1) 

for some Z^r vector n a , rank of the gauge group of the 2d GLSM is given by 

K = ^(n a -l) = ^n a , (4.2) 

a=l a=l 

masses of the 4d theory and twisted masses of the 2d theory are related to each other in the 
following way 

M a = m a - |e , M a = m a + |e , (4.3) 



3 We change the notations of 11 to make them consistent with our notations. 
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and coupling constants obey 

t = t + \{N + 1). (4.4) 
Quantitatively the CDHL duality states that in NS limit, the chiral rings of the 4d and 2d 



theories are isomorphic and we can relate effective twisted superpotential (1.1) from the 4d 



gauge theory with the effective twisted superpotential of the corresponding 2d GLSM as 
follows 

W(0 Q =m a - n a e) - W(0 a = m a - e) = W*£(n a ) . (4.5) 



This means that the sets of stationary points (vacua) of the two superpotentials are isomor- 



phic and the above equality holds in the corresponding vacua. The equality (4.5) has been 



proven by computing the Nekrasov partition function on the Higgs branch of the theory (4.1 ) 



and deriving the effective twisted superpotential then matching it with the 2d superpotential 



on shell 11,28 . Let us mention a related contribution 29 , where moduli space of vortices 



was shown to be a submanifold of the instanton moduli space of the Omega deformed 4d 
theory. 

Below we shall provide some further supports for this duality based on the study of non- 
Abelian BPS vortices in the Omega deformed four dimensional theory. This will be done by 
identifying the classical 2d theory living on the vortex along the course of the Shifman-Yung 
program [To]. 



4.2 Constructing non-Abelian vortices Nf = N 

Let us now construct the action of the M = 2 SQCD in the NS Omega background in four 
dimensions. Although in 11,28 the superconformal Nf = 2N case was considered, from the 



viewpoint of the non-Abelian vortices it is more instructive to start with the left boundary 
of the stability window Nf = N. 



Action. Let us begin again with the undeformed SQCD Lagrangian in four dimensions 
with N = N f 



C 



Air 



-Tr 



ra\2 



+ I d^Oi Q ie y Q l + Q i e~ v Q i ) + I d 2 6( Q&Q* + m^Qi + H.c 



(4.6) 



where V is an adjoint Af = 1 vector superfield, $ = ((f), A, D) is an adjoint SU(N) C U(N) 
chiral superfield, quark superfields Qi = ipi, i^) and Qi = {q^ijj^Fi) are transformed in 
N and N of the SU(N) and as Nf and Nf under global flavor SU(Nf) respectively; rrii are 
quark masses and r = ^ + ^ is coupling constant. 



Omega Deformation. Our task now is to construct the Omega deformed theory. As 
in the pure SYM case the deformation can naturally be understood in terms of the six 
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dimensional N = 1 theory 16 . It is convenient to use dual frame description Gab = ^a^b ■ 
The components of sixbeins read 

e ( m)=C; c M = n?, C W = 0, e^ = S a b . (4.7) 
Using the above equation we can rewrite the kinetic term for squarks 

e { a B) V {B) q = V a q-in™V m q. (4.8) 

Thus we have 



|VaC = |V m gf + 1(0 - ^ m V m )gr , (4-9) 

and analogously the kinetic term for anti squarks. The bosonic part of the action after quark 
masses m; and rhi are included reads 

£ = w F L + ^iv m - F mn n n \ 2 + ^1 <K0 - iv m (n m r - n m r) + 9 2 (qr a q - qr a q)\ 2 
+ ||V m g| 2 + ||V m g| 2 + §|(0 - - i^V^ftl 2 + ||(0 - m, - ^ m V m )g,| 2 
+ 2g 2 \qr a q\ 2 + - N£ FI \ 2 + £ (|g| 2 - |g| 2 ) 2 , (4.10) 

where we have included Fayet-Iliopoulos term This theory has U(N) C x SU(N)f global 
color and flavor symmetry group. 

Supersymmetry transformations. In what follows it is convenient to package squarks 
and anti-squarks into a single vector tp = (q l , q 1 ), where / = 1 for squarks and / = 2 for 
antisquarks. Supersymmetry acts on the fields in the following wajj^] 

6<j> a = -V2( af (K f - ^K f ) + <& a6l \% , 

5r a = -iV2CfV a aq if + 2iCl [Ur a ))q J f ~ <^ 4 V^gj) , (4.11) 

where z = 1, . . . , N runs through fundamental representation, a = 1, . . . , N 2 runs through 
the adjoint representation of U(N), f,g = 1,2 denote SU(2) R-symmetry index, and the 
D-term contribution in the first line above has the following form 

D»' = -0r>^ - g 2 (q g r a qf - E a /) , (4.12) 

where the generalized Fl-term reads 

S*' = ^V^O^r - Sr*P)8' a + & Ig 5 a N2 . (4.13) 

The first term we have already seen in the previous section, similar story here - it is generated 
by the Omega background. The second contribution to 5 is the standard FI term. Here, 
we formally kept the FI parameter £ f FI triplet. Usually only diagonal part 



as m 



30 



We use chiral notation here. 
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is left over, as it simplifies the calculations, however, it is absolutely unnecessary. As we 



can see, generalized FI parameter (4.13) is a sum of the two terms, the field dependent FI 



term, which appears due to Omega deformation, and the conventional U(l) FI term, which 



is normally considered in super symmetric QCD. As we know 31 , presence of the latter does 



not affect the supersymmetry of the theory, however, it's broken to the (2, 2) SUSY due to 
the former. Similar to the pure SYM case, considered in the previous section, generators 



(2.14) form the supersymmetry algebra. We shall use this information to extract the BPS 



equations later. 

Classical vacua. Vacua of the theory can be chosen similarly to the undeformed SQCD. 
Indeed, after identifying 

(j) a = m a q ai = q ai = y/^S^ , (4.14) 





we find that the potential in (4.10) vanishes. Alternatively one could have put q l 



and work only with squarks q l (see 10 for details). This vacuum is invariant under the 
color-flavor rotations 



G c+f : cP^U^Uf, ql^(U- l ) b a q{U f 



3 J 



(4.15) 



where U c = Uf G G c+ f. In the most generic case, when masses m a are arbitrary, one has 



G c+f = S{U{n x ) x ■ • ■ x U(n k )) , £ 



n ; 



N. 



(4.16) 



i=i 



where S(. . . ) stands for the stabilizer. If all masses are different then G c+ f = S (U(1) N ) = 
rj^g p a ttern of the symmetry breaking depends on the relationship between the 
masses rrii and the FI parameter £pj. In what follows we shall assume £ 3> mf for any i and 
all the masses to be of the same order, thus in our case we have the following breaking 



U(N) C x SU(N) f ^ U(l) x SU(N) c+f ^ U(l) x G c+f . 



(4.17) 



Note that U(l) factor in the above formula will be very important in constructing vortices, 
which we shall now do. A vortex configuration will further break the above symmetry 
in a nontrivial way - the above mentioned U(l) will be coupled to the generators of the 
Cartan subalgebra of G c+ f. By using the extra U(l) symmetry we can also change the 



second equality relation in (4.14), and this is what exactly is done in the DHL paper. For 



the convenience of the calculations of 11 squarks have charge —3/2 and anti-squarks have 



charge +1/2 with respect to this symmetry. In the current paper it is more convenient to 



keep the condition (4.14) on the vortex solution as well. 



Vortex configuration. In order to find non-Abelian BPS strings we need to organize 
winding around the z-axis. Thus we allow one of the flavors, say q N , to depend on the 
azimuthal angle e mip q , where n is an integer. Algebraically it corresponds to breaking the 
symmetry of (4.17) down to C/(l)diag x SU(N — 1), where the first £/(l)dia g factor is the 
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diagonal subalgebra of the U(l) from (4.17) and the N — l'st Cartan generator of the G c+ f. 
Then for the two terms in the second line of (4.10) read 



VD ^-m^ + hed'+epiA^q^] 2 . (4.18) 

It vanishes provided that the expression in the parentheses above is equal to zero. So we put 

<f>% = m N -ne-ep{A v )%. (4.19) 

It can certainly be generalized to the case where more squark fields have angular dependences 

(f) a = m a - n a e - iQ m A a m , (4.20) 

where n a is an integer valued vector, which is intended to count the number of flux quanta 
which flow through the vortex. We see that the above classical vacuum equation related the 
adjoint scalar and the gauge field. 



Vortex BPS equations. While studying a 1/2-BPS object we work with the half of 
supersymmetry algebra which acts trivially on it. In the case at hand this algebra is generated 



by (2.14). Remarkably it coincides with the BPS subalgebra of the non-Abelian vortex 
considered by Shifman and Yung 10 . Thus even in the Omega deformed background in the 



NS limit, the vortex configuration we are considering in this section will remain 1/2-BPS. 



Performing Bogomol'ny completion of the action (4.10) we get the following energy den- 



sity 



C 



i 

2g 2 



\(B a 3 ) 2 + g 2 (qr a q-E a ) 



+ ^|(V 1 + *V 2 )r-(^2-^i)^| 
|(Vi + iV 2 )q\ 2 + N£ FI B^ + ±d m (n r > " 



v " a -n m (j) a )Bi . 



(4.21; 



Here we assumed that the adjoint scalar and gauge field are only aligned along the Cartan 
subalgebra of the gauge Lie algebra. The last two terms in the second line of the above 
expression are total derivatives, but due to a different reason: the former is the Abelian field 
strength, which gives circulation of the gauge field after removing one integration, the latter 



involves cL derivative and is of the same kind as (2.16). We can see that the Lagrangian 



(4.10) under the constraint (4.20) and color-flavor locked condition q 1 = q 1 takes almost 



exactly the same form as for the undeformed case considered by Shifman and Yung 10 . It 
means that the BPS construction for the vortex will also be almost exactly the same. The 
only difference is that adjoint scalar a will have a nontrivial profile defined by the magnetic 
field and the Omega background. The corresponding BPS equations read 



5 3 a + 5 2 ferY 



0, 



(Vi + iV 2 



(Vi + iV 2 y = 
(0 2 - iVtx)Bl = 



(4.22) 



where, again as in (4.13), the color index a = 1, . . . ,N 2 runs through all U(N) generators. 
For convenience we can split up U(l) and SU(N) parts and rewrite the first equation above 



using the definition of the generalized FI term (4.13) 

9 2 (\Q\ 2 



B* 



Bl + g 2 (qiT a q*-C F l) 



0, 
0. 



(4.23) 
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where we have denoted 



aa.a 



tt aa (h a ) 



ech a + ec 



,N, 



(4.24) 



the non-Abeilan FI field. In its absence equations (4.23 ) and second equation in ( 4.22[ ) exactly 
reproduce the BPS set considered in 



10 



hence the solutions for the profile functions can 
be extracted from there directly. Thus the modification to the BPS vortex equations in the 
Omega background consist of introducing £%fi (which for some configurations can vanish) 



and the nontrivial profile for the adjoint scalar dictated by the third equation in (4.22) 



Asymptotic behavior of solutions. Let us for the moment assume that <p a is invariant 
under rotations around the z-axis, in other words £® FI vanishes. From the analysis of the 
previous section we conclude that it happens when does not depend on the azimuthal angle 
ip. Then, we know the solution for the magnetic field in all color directions, since it is exactly 
the same as in |10|. In particular, far away from the vortex, the gauge field exhibits 1/p 



behavior. This makes the quantization condition (4.20) physical and well defined. Indeed, 



it tells us that the adjoint scalar at large p approaches its vacuum value 



m 



e(n a + k a ) 



(4.25) 



where k a is a Ztv- valued vector of winding numbers of along the different Cartan color 
directions. We still need to figure out what k a is in terms of n a . The reasoning for that comes 
from the following physical requirement - string tension (energy per unit length) should be 
finite. Indeed as in the undeformed case, the conclusion comes from the requirement that 
|V m g| 2 terms are finite. 



+oo 



dpp\V m q 



1 1 2 



(4.26) 



This was archived by a proper asymptotic behavior of the azimuthal component of gauge field 
such that the integrand above could decay fast enough. Indeed, let's say q ~ e mip q(p), 
thus the integral becomes 



dp - I (in — iA^p) q 
P 



i\2 



(4.27) 



so Atp — > n/p at large p. In other words, A v should be proportional to the number of flux 
quanta which flow through the vortex. So, given (4.25 ) we can easily figure out that k a = —n a 
and tends to its undeformed value m a at large radial distances. Thus we conclude that 
the adjoint scalar interpolates between 



at p = and 



at p = oo, see Fig. [3j 



h a = m a - n a e 



a = m a 



(4.28) 
(4.29) 
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Figure 3: Scalar field interpolating between two different minimum values of the potential 



Now we understand why in the left hand side of (4.5) involves the difference of the 
superpotential in two points - they are two different minima of the effective potential for 
in four dimensional theory, and the vortex can now be viewed as a kink which interpolates 
between these minima! Evidently (4.5) is only applicable for nonzero e. Note that the 
superpotential in the left hand side of (4.5) is evaluated at values of <fi which are shifted by 
a unit of e. As we shall explain below it happens because of an additional U(l) twist of 
(anti)squark fields. 

Let us mention that a completely different situation occurs if (p acquires nontrivial topol- 



ogy in the spirit of the previous section. Then BPS equations (4.22) do not decouple any 



longer. We expect a significant change in the asymptotic behavior of the solutions in that 
case. Investigations in this direction will be reported elsewhere. 



4.3 Constructing non-Abelian vortices N < Nf < 2N 



Let us now address semilocal vortices [32]. In order to be more generic we shall keep N = 
Nf — iV generic inside the conformal window. The vacuum condition is generalized as follows 



m . 



Qa 



S&j8< a , Ki<N 
0, i > N > 0. 



(4.30) 



Symmetry breaking pattern is similar to (4.17), only there is a residual global symmetry left 
due to the additional quark fields 



G c+f = SU(N) c+f x SU(N) . (4.31) 
Most recent review of the semilocal vortex constructions can be found in [33] . The structure 



of the BPS equations (4.22) will not change, only the flavor index will now range i = 
1, . . . , N + N. On the level of the low energy effective action the vortex theory will be 
modified by adding iV so-called size moduli. Kinetic terms of size moduli bring logarithmic 
divergence to the energy density of the effective theory, therefore one has to introduce an 
infrared cutoff. 
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GLSM description. As usual, a GLSM description of the theory is more effective for 
computations. The 2d theory which is dual to the 4d SQCD in the NS Ome ga b a ckgr ound 
with Nf = N + N quarks is given by the following Lagrangian provided that (4.1)-(4.4) hold 



£ = Tr / d A 6 



2e 2 



+ Tr / d 2 6ri: + H.c. : 



1 N N 



i=i 



i=i 



(4.32) 



where the trace is taken over the adjoint representation of U(K) gauge group, <3> is adjoint 
chiral multiplet, and £ is field strength for 2d vector superfield V. The second line in the 
above Lagrangian represents the twisted F-terms of the theory. There are N + N +1 twisted 
mass parameters turned on including N + N masses for X and Y fields together with the 



twisted mass for the adjoint scalar <3>, which according to 11,28 equals to e. In the limit 
e — > oo the gauge field becomes non dynamical, and we can integrate it out. In this limit 
we can recover the geometry of the NLSM's target space, which naturally appears in the 
derivation of the low energy theory. 



In order to get the effective twisted superpotential in the right hand side of (4.5) we need 



to integrate out X's, K's and $'s in (4.32). When Nf = 2N C the theory is superconformal, 
the coupling does not run and no dynamical scale is generated. 



~ K N (\ M\ K N (\ 

C(A)= e EE/ ^ -EE/ - 

a=l i=l ^ ' a=l i=l \ 

\ K 

' o=l 



A„ - Mi 



A a — Ab — e 



(4.33) 



where f(x) = x(\ogx — l). Note the change of the coupling constant to f compared to (4.32). 
Minimizing the above superpotential we arrive to the ground state equations 



N 



n 



A . ; - M t 



K 



i \j - Mi 



TT Aj - A fc - e 



(4.34) 



which coincide with Bethe ansatz equations for the twisted anisotropic Heisenberg SX(2,R) 
magnet. This observation quantifies the so-called Bethe/gauge correspondence for the Af = 2 
SQCD. 

Theories with N < N can be obtained from the conformal theory by sending some masses 
to infinity and renormalizing the coupling constant. Dynamically generated scale Aqcd will 
then appear. 



A note on the Af = 2* theory. Recently a similar to DHL and CDHL study of the 
M = 2* theory in five dimensions appeared in the literature [34], and a duality with a 
three dimensional integrable system was discussed. Although the calculations involving the 
Nekrasov partition function look very similar to 11,28 , there is a technical difference: the 
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Higgs branch condition of the 5d theory looks similar to (4.1), however, there is no shift by 
N in the rank of the 3d gauge group like in (4.2). Clearly, is occurs because fundamental 
matter in 1128 and adjoint matter in 34 contribute differently to the Nekrasov partition 



function; nevertheless physical understanding of the second duality remains to be uncovered. 
Since there are no BPS vortices in M = 2* theory, one cannot apply the method we used in 
the current section. 



5 Brane Constructions and Dualities in Integrable Sys- 
tems 

Solutions to the Bethe ansatz equations mentioned above correspond to the ground states 
in the world volume theory on the non-abelian strings realized as D2 branes, hence it is 
desirable to translate the full powerful machinery of the integrable systems into the brane 
language. In this section we focus on the particular issue namely the realization of known 
dualities between quantum integrable systems using brane language. 

We shall first review the Hanany-Witten type IIA brane construction which yields the 
M = 2 SQCD and integrable systems related to it - the XXX spin chain and Gaudin 
model together with the dualities these models are involved in. Employing the Gaudin/XXX 
duality we will be able to give a vortex interpretation of the AGT duality in the next section, 
where the XXX model appears on the M = 2 theory side and the Gaudin model naturally 
arises in study of Liouville CFT. Here we shall make some preparations to that study. 
In addition to that the Gaudin/XXX duality will be examined by studying the simplest 
examples of Argyres-Douglas type points and wall crossing phenomena in presence of the 
Omega background. At the end we shall discuss yet another duality between spin chains 
and Calogero-Moser systems. 



5.1 Dualities from the Hanany-Witten brane construction 

Brane configuration for the M = 2 SQCD employs the Hanany-Witten construction |35|. As 



it was shown in 11 and further explained here in Sec. 41 in presence of Omega background 



the Higgs branch condition gets deformed (4.1 ). Hence the positions of the flavor D4 branes 



are shifted by n a e for each color (see bottom picture in Fig. |4J). It contains two NS5, N D4 
branes which are stretched between the two NS5 branes and two sets of semi infinite D4's 
which are attached to NS5's. All D4 branes occupy 01236 directions, NS5's lie in 012345 
directions. 
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Under geometric transition the brane configuration described in 11,28 interpolate between 
the 4d theory and the 2d theory. The latter can be obtained by moving the right NS5 brane 
in the 7th direction and emerging D2 branes (037) which are stretched between this NS5 



23 







m\ — n\€ 




mi 


m\ 

m 2 - n 2 e 
m 2 


rh 2 






m N 




jtin — n^e 


m N 







D2 






fn 


L 






/ h 


2 


rh 2 




/ 








A'N 










m N 



Figure 4: Top: Type IIA brane picture. Positions of semi infinite D4 branes in 45 plane is given 
by rrii + n^e, where i = 1, . . . , N. In the deformed configuration the Higgs phase of the theory is 
given by the condition cti = rrii + rue. Bottom: 2¥ortex strings as D2 branes stretched in the 7th 
direction. 



and D4's (see right picture in Fig. |4J). The value of xj gives tension of D2 strings which is 
equal to e in our construction. 

The rank of the gauge group of the two dimensional GLSM is given by summing up all 
the D2 branes K = J^rU, where, we remind, hi = rij — 1. The low energy dynamics of 
the two dimensional theory is given by the effective twisted superpotential and the following 
ground state equations 

^j^wr'U^^' (5 ' 1) 

which is the Bethe ansatz equations for the anisotropic SL(2) spin chain. Note that for 
generic 2d masses M a and M a at each spins at each site a = 1, . . . , N have different repre- 



sentations. Indeed, in order to match each term in the left hand side of (5.1) with phases of 
anisotropic chain 

where v a are anisotropies and S a are spin^j one identifies [5] 

M a = 6 a -S a e, M a = 6 a + S a e. (5.3) 

5.2 The Gaudin/XXX duality 

It is known that the Gaudin model [36 enjoys several dualitie^] First we recall the duality 
introduced at the classical level in [37 . It relates the rational Gaudin model with SL(N) 
group at M sites and SL(M) group at N sites. The positions of marked points Zi on the 
sphere corresponding to the inhomogenities and the diagonal element of the twist matrix 
get interchanged. At the classical level the spectral curves and the action differentials are 
equivalent. At the quantum level the Bethe ansatz equations reflect this symmetry at the 
level of spectra. 

Let us explain this symmetry in the brane picture. Let us first remind ourselves the 
similar symmetry in the Toda system discussed in j3j. It the Toda case this symmetry 
merely implies the equivalence of 2 x 2 and N x N Lax operator representations which can 
be explained as the 90 degrees rotation of the viewpoint of the brane picture. In the first 
representation the gauge group is connected to NS5 branes, while in the second case it is 
defined by the number of D4 branes in the IIA picture. 

If we add the fundamental matter and consider the conformal case there are additional 
data which have to be matched via the duality. In the 2x2 representation the SL(2) twist 
matrix emerges which reflects the positions of NS5 branes in the 6-10 plane Fig. [5] The 
masses of the fundamentals provide the inhomogenities at the corresponding lattice sites. 
Upon the 90 degrees rotation similar to the Toda case the two sets of data get interchanged. 

The duality between a pair of rational Gaudin models can be generalized to a similar 
duality between a trigonometric Gaudin model and a XXX spin chain via the so-called 



5 For example, S a — — 1/2 gives the SL(2) chain 

6 Some details about the Gaudin model are given in App. [a| 
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Figure 5: (6 + HO, 7) section of the HW brane construction (view from "below"). 



qI(M)/qI(N) duality 38 . For M = N = 2 Bethe ansatz equations read as follows^] 



Mi - M 2 - e + y> u b e y-^ 2e 



for trigonometric Gaudin, and 



6=1 



«2 



(5.4) 



n 



Z-2 



"2 



TT Aj - Aj - e 



1 Ai + A^ a + « a e 1 Ai - A,- + e 

a=l j=l J 



^2 



(5.5) 



for the SL(2) XXX chain. The Mukhin-Tarasov-Varchenko (MTV) duality 38 states that 
(5.4) as set of equations with respect to ti, . . . t K2 has isomorphic space of orbits of solutions 
with the one of (5.5) as set w.r.t. Ai, . . . , \ U2 provided that 

k 2 + k 2 = v x + v 2 . (5.6) 



Parameters M.\ i2 and z X2 are generic. We can now recognize (5.1) in (5.5) with 

M a = -M a , M a = -M a -K a e } K = v 2 , N = 2, z x = l, z 2 = q , (5.7) 

and parameters k 12 and V\ will be specified later. Also it will be more useful for us to use the 
4d masses instead of the 2d ones. We can then rewrite set of MTV dual equations ( 5.4|5.5 ) 
as follows 



-mi + m 2 — e v x e 



Ke 



t{ — Z\ ti — z 2 



n 



a=l 



2e 

j=l ~~ 3 
K 



Aj - m a + |e _ z 2 tt Ai - Xj - e 
\i-m a -\e z x \\\i- \ j + e 



2 e z\ j=1 ,i, 



(5i 



7 We have adopted the notation and made some change of variable compared to 



38 



26 



Thus we can see that twists z±, z 2 , corresponding to the positions of the NS5 branes in 6-10 
plane Fig. |5j and masses of the fundamentals mi, m 2 interchange their roles upon the duality. 
We see that matching to the BAE corresponding to U(2), Nf = 4 SQCD shows that the 
strange nonequal mass shifts to the fundamentals and antifundamentals (4.3) have now clear 
interpretation within the duality. Namely, the number of the Gaudin Bethe roots yields the 
asymmetry between the fundamental and antifundamental masses. Also Gaudin spins match 
with the number of Bethe roots at the XXX side. Later in the next section we shall use 
these spins in order to make the AGT duality manifest. 

Let us emphasize that the Hamiltonian of the Gaudin model is nothing but the r.h.s. 



of the Knizhnik-Zamolodchikov (KZ) equation 39 on the sphere with L + 3 marked points 
* |40) 

b 2d ^=U Gaud y( Zi ), i = l,...,L, (5.9) 

where b is some constant. In the next section, when we will discuss Liouville theory on the 
same Riemann surface, we shall specify its value. 

One could also introduce the so called dynamical operators with respect to boundary 
conditions [38]. Under the bispectral duality transformations the Gaudin KZ operator and 
the dynamical operators get interchanged as well. The number of marked points in the 
N x N representation of the Lax operator corresponds to the number of NS5 branes involved 
in the gauge theory brane construction. 



5.3 Bispectral duality and Argyres-Douglas points 

Classically, the bispectral duality just states that two systems have almost the same (differ 
by the simple factor) spectral curves and action differentials. At the quantum level the 
situation is more subtle. Since naively the bispectral duality connects the systems with 
different degrees of freedom, one should be able to analyze the phenomena of merging of two 
degrees of freedom into the single one. Below the simplest example shall be considered. 



Let us now consider (5.4) again, this time we identify Aii = —M. 2 = le- Then one has 



K-2 



21 — 1 V\ V 2 v — • z, 

+ TZT7. + TZT. ~ Yl T~T. = °> • = !.-.«.. (5.10) 



ti ti Z\ ti Z2 . ^ ti tj 



for the Gaudin model and 



l-e s a + l-e z 1 ~p l s a -s b -e 

I — = 1, a=l,...,V2 (5.11, 

J- J- e — «, 4- e 



s a - I - e - Kie s a + I - e - k 2 z 2 * s a - s b 

0=1 

for the XXX model. Integers K a , v a satisfy the relation K\ + K 2 — V\ + v 2 . One of the results 
in 



38 is the precise correspondence of the orbits of solutions to the Bethe equations under 
the group of permutations of variables (permutations of ti, ...,t K2 for the Gaudin model and 
of S\, s„ 2 for the XXX model). At first glance such correspondence seems to be quite weak 
as it does not establish a direct connection between the roots of both systems and does not 
allow one to simplify one set of equations having known the solution of the other. But it 
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preserves one important feature of the XXX model, namely the degeneration locus which 
could be called a little bit loosely "quantum Argyres- Douglas (AD) points' : 



41 



The "classical" Argyres-Douglas manifold is the locus in the moduli space of the theory 
where different vacua merge together. To get such merging typically one starts with M = 2 
theory and perturbs it down to M = 1 and tunes the external parameters like couplings and 
masses. In the "quantum case" we have one more parameter from Omega deformation e 
and the AD manifold involves this additional coordinate in the parameter space. The rest 
is the same and AD manifold corresponds to the colliding vacua. In this subsection we shall 
normalize e—1. 

Since the solution to the BA equation correspond to the vacuum state "quantum AD 
point" corresponds to the appearance of multiple roots. We shall consider a simple case 
fti — v 2 = 2, k 2 = v\ = 1 as an example and find the AD manifolds for both models. Bethe 



ansatz equations the XXX chain (5.11) then read 



si 



I — 1 Si + I — 1 Si — s 2 — 1 



«1 

s 2 



I -38! + I 

l-ls 2 + l 



2 Sl 
ls 2 



S2 
Si 



s 2 - I - 3 s 2 + I - 2 s 2 - si + 1 



q, 



(5.12) 



and the Gaudin system is described by a single Bethe equation. In general, the XXX BAE 
contain a number of degenerate solutions, such that for some i,j the roots coincide Si = Sj. 



The vacua of the theory correspond to the non-degenerate solutions of the system 11 



Obviously every solution to the degenerate system is a solution to the full BAE system. 
This property can be used to lower the degree of the Bethe equations. In our case the 
degenerate solution is si = s 2 , so from (5.12) we obtain 



s-l-ls+l-1 

s -I -3s + 1-2 



-q- 



(5.13) 



We can now solving the first equation of (5.12) with respect to s 2 and substitute this solution 



into the second equation to obtain a polynomial roots of which solve the XXX BAE system. 



In order to eliminate the degenerate roots we merely need to divide this polynomial by (5.13 ) 
To find the AD manifold we calculate the discriminant of the reduced polynomial 



D x (l, q) = 4(1 - qfq 2 (Al 2 q 2 - 8l 2 q + 4l 2 - Alq + 41 + 8q + l) 2 
(4/V - 32/V + 56/V - 32/ 2 g + Al 2 + Alq 4 - 3Qlq 3 
+ 28lq 2 + Alq + q 4 - 18g 3 + 17q 2 - 8q) . 



(5.14) 



However, the set D x = still contains extra roots. Although the equation was divided by the 
degenerate one, the discriminant still captures the cases when the roots of the reduced system 
coincide with the roots of the degenerate system. In order to exclude such cases discriminant 



D\ (5.14) must be divided by the resultant of the reduced and degenerate polynomials. This 
resultant turns out to be precisely the polynomial in the last parentheses in (5.14). Thus 



the resulting AD set is described by the zero locus of the following polynomial (we do not 
discuss trivial cases when q = and q = 1) 



2„2 



D(l,q)=4l 2 q 



8l 2 q + Al 2 



Alq + Al + 8q + l 



(5.15) 



28 



which is precisely the discriminant of the Gaudin equation written in a polynomial form! 

The above example describes the way to calculate the quantum AD set for the XXX 
BAE. First, we identify the degenerate subset and divide all equations by the corresponding 
polynomials. Then we find the discriminant of the reduced polynomial and divide it by all 
possible resultants with the degenerate polynomials. The same procedure can be done for the 
Gaudin system. The resulting polynomial describes the set of the quantum AD points and 
coincides for XXX and Gaudin systems. Certainly we have considered the simplest example 
and this issue deserves a separate study. 



5.4 Walls of marginal stability 



One more issue we shall briefly discuss concerns the interpretation of the walls of marginal 
stability in the Omega deformed theory in the language of the quantum integrable system. 
The wall can be described in terms of the superpotential however the Yang- Yang function in 
the quantum integrable system has the interpretation of the twisted superpotential as well. 
Hence we could formulate the problem of finding the walls of marginal stability in terms of 
the YY function for the spin chain or Bethe ansatz equations. 

The e-deformed theory has exact effective twisted superpotential (4.33). XXX Bethe 



ansatz equations (4.34) specify the positions of the vacua. Kinks in the two-dimensional 



theory which interpolate betw een these vacua possess two kinds of charges: Noether charges 
Mi, Mi or 2d twisted masses (5.3) and the topological charges which are given by the dif- 



ference of the vacuum values of the superpotential (4.33) evaluated at the vacua the kink is 
hopping between. A wall of marginal stability is the solution to the following kinematical 
condition of the decay of the kink into its constituents 



3m 



^top 



^Nothor 



Mi - M 2 



0. 



(5.16) 



We follow the procedure described in 42 . The idea of the analysis is that the vacua and 



the vacuum values of the superpotential should be continuous functions of all parameters 
across the wall of marginal stability. We investigate the system depending on a single complex 
parameter 9 and make all the deformation parameters to be equally spaced on a circle 



6i = 9 exp 



2nil 



(5.17) 



where 9\ are the spin chain impurities and contribute to 2d masses (5.3). Note that the 



above Z^ choice of 0's is a simplification as it certainly cannot be done for generic masses. 
However, as it argued in 42 , it is generic enough for the study of wall crossing phenomena 
in 2d. 



not as sep- 



We can now interpret the vacuum solutions to the BAE equations 

arate functions, but rather as branches of some continuous function \ vac {9) on the complex 
plane. The number of branches £ is of course finite and coincide with the degree of the BAE 
system as a system of polynomial equations. The superpotential is the continuous function 
depending on \ vac {9). It has infinitely many branches but the sequence of branches has a 
certain periodicity of order I. 
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The supersymmetric 1 theory considered in 42l 43 has infinitely many walls of 



marginal stability The vacuum as the function of the twisted mass parameter Mq 

M, = M exp(^^ (5.18) 

has N branches and the branches of the superpotential differ by the quantity proportional 
to m 

Wffl(M ) = exp (VV(°2(M ) - 2vrmM ) . (5.19) 

The condition for the kink interpolating between nth and (n + l)st vacua to merge with a 
state of Nother charge mi — mo on the wall of marginal stability reads 

[Mac 2mnexp(2g)\ 

fHe — /9 , N 1 = 0. 5.20 

^ Mi exp V ' 

The above equation describes N walls which are roughly logarithmic spirals with infinitely 
many branches. 

In the case of the Omega deformed theory the situation changes a little bit. In order 
to depict the walls of marginal stability graphically in general case, one needs to solve the 



system of BAE. We take the case iV = 3, K = 1 in (5.1 ) as the simplest example. The BAE 
system is 

(A — St) 3 — 8 3 ~ [b ' ZL> 

The solution has three branches. The difference between nth and (n + 3)rd branches of the 
superpotential equals 4ni. The walls of marginal stability are depicted in Fig. [6} In the limit 
e — > the XXX system passes into the Gaudin system. The walls of marginal stability can 
be plotted for the Gaudin superpotential too. They possess the main features of the walls 



of marginal stability for the XXX system. The superpotential (4.33) at small e becomes 

K K N N 

W G audin(A,) = G X j ~ 2 Sl lo S ( X i ~ 9 ^-J2 lo S ( A i - ' ( 5 - 22 ) 

3=1 j=l 1=1 j^k 

where G is external field. In the case iV = 3, K = 1 the solution to the corresponding 
Gaudin equation has three branches. The difference between nth and (n + 3)rd branches of 
the superpotential remains Ani as in the XXX case. The corresponding walls of marginal 
stability are drawn in Fig. [7j 

5.5 On the spin chain/ Calogero duality 

One may wonder if the bispectral duality between two different types of the spin chains 
has any relation to a similar duality discussed in the context of the Calogero-Ruijsenaars 
family of the integrable systems. The answer turns out to be positive. The prototype of this 



duality has been discovered in 44 , where the correspondence between the zero locus of the 



Toda Hamiltonians in the phase space and the manifold associated with the quantum model 
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Figure 6: The walls of marginal stability for the XXX superpotential on the 6 3 plane. The param- 



eters are: q 



1, S = |, e = 1. The AD point is at 9 3 



has been discovered. This correspondence has been generalized for the rational Calogero 



model 45 at fixed coupling which turns out to be dual in the same sense to the rational 



Gaudin model |46|. 

The identification of the parameters goes as follows. The inhomogenities in the Gaudin 
model Zi are identified with the coordinates of Calogero model Xi, while the eigenvalues of the 
Gaudin Hamiltonians Hi are identified with the Calogero momenta. The inhomogenities are 
Poisson dual to the positions of the marked points, hence their identification with the Poisson 
pair in the Calogero model is natural. The duality above can be extended from the zero locus 
at the Calogero side to the arbitrary values of the classical Calogero Hamiltonians |47| 



n 



Cal,k 



0*(A) 



(5.23) 



where cr& is the k-th symmetric power of the Lax connection eigenvalues. It turns out that 
the Lax eigenvalues at the Calogero side are mapped onto the eigenvalues of the twist matrix 
at the spin chain side. 

As we have discussed above the rational Gaudin model enjoys the marked points/twist 
duality and in some sense is selfdual. Its bispectral dual - the rational Calogero model 
is selfdual as well in the same sense. The bispectral duality can be generalized to the 
trigonometric and relativistic cases 48 50 . Thus the trigonometric Calogero-Moser model 
is known to be dual to the rational Ruijsenaars-Schneider model |51,49^J The quantum 



'See also 52 
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version of this duality has been elaborated in 53 . Recall that in Sec. 5.2 we discussed another 



bispectrally dual pair between the trigonometric Gaudin and the XXX models [54]. We can 
now see that the bispectrality at the Calogero-Ruijsenaars side matches perfectly with the 
duality at the Gaudin-XXX side. Note that this duality has the clear interpretation in terms 



of the Chern-Simons theory with inserted Wilson lines and its Yang-Mills degenerations 55 
The relationship between the two dualities is summarized in Fig. [8] 

It would be interesting to make the next step and consider the selfdual trigonometric 
Ruijsenaars model at the Calogero side of the correspondence. Its dual on the spin chain 
side is expected to be the XXZ chain which, according to the correspondence considered 
above, should enjoy some kind of bispectral selfduality. Another important issue concerns 
the generalization of these dualities to the elliptic integrable models. Not much is known 
about the self-dual elliptic model yet 



55,56,12,57 



(see, however 

To conclude this section note that curiously we arrive to the claim that since the twists at 
the spin side corresponds to the spectrum of the dual Calogero model, therefore the SQCD 
gauge coupling, which corresponds to the twist, could be quantized! This point deserves the 
further analysis. 



6 The AGT Correspondence in the NS Limit 

The Alday-Gaiotto-Tachikawa duality [15] relates a conformal block of the Liouville CFT 
on a Riemann surface of genus g and n punctures with the instanton part of the Nekrasov 



32 



Trigonometric 
Calogero-Moser <Cj 
model 



A 



bispectral 
duality 



Rational 
Ruijsenaars- 
Schneider model 



MTV 



V 

Twisted anisotropic 
XXX chain 



bispectral 
duality 



Trigonometric 
Gaudin model 



Figure 8: A pair of bispectral dualities mapped onto each other 



partition function of a quiver gauge theory naturally associated with this Riemann surface. 
Furthermore it concludes that the n-point correlator in the Liouville theory on this Riemann 
surface can be evaluated as an integral of the square of the absolute value of the full Nekrasov 
partition function. 

Having done the above analysis on spectral duality between integrable systems, we arrive 
at an interesting observation which envisages the AGT correspondence for the Liouville 
theory with large central charge on S 2 with four punctures and U(2) SQCD with four flavors. 
It has already been addressed in the literature earlier [58] where a Liouville conformal block 
at large c simply becomes a hypergeometric function of the conformal dimensions and the 
instanton number (59) was used. Then the authors figured out that only chiral terms in the 
Nekrasov partition function will contribute, either e\ or €2 have to vanishes; it enabled them to 
identify each multi-instanton contribution with corresponding terms of the hypergeometric 
function's expansion. The proof is rather formal and it will be desirable to have a more 
physical rationale to it. The current section is intended to fulfill this goal. 

In order to see how the AGT relation comes about, in what follows we shall relate both 
the gauge theory and the Liouville theory to a pair of integrable systems which enjoy a 
certain duality between them. The roadmap we shall use to guide us through this section is 
presented in Fig. [9j Starting from the 4d gauge theory on the top right of Fig. [9] we shall 



use the results of Sec. |4jand 28 in order to relate the 4d theory with the corresponding 2d 
GLSM. As it was shown in |28|, the equivalence was established to all orders in the instanton 
parameter. Thus, instead of working with each summand of the instanton partition function, 



as all known proofs of the AGT 58 



60 did so far, we shall look at the entire expression and 



the effective twisted superpotential (1.1 ) which follows from it. Following this line of thought 



we provide a physical rationale for the AGT correspondence in the NS limit, namely, modulo 
certain dualities between two given integrable systems, it is reduced to the 4d/2d duality 
in the NS limit |28|. The latter exists due to non-Abelian semilocal BPS vortices which we 
have discussed in Sec. |4j also vortices exist provided that the corresponding FI parameter is 
turned on, hence we uncover why an extra U(l) factor on the gauge theory side is important. 
The quantum duality, together with its brane interpretation, was discussed earlier in 



Sec. 5.2 Now we shall start with the left column of Fig. M by reminding ourselves how 



the Gaudin model is related to Liouville conformal blocks, and later on, by means of the 
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Figure 9: Roadmap of the AGT duality in the NS limit. The main statement (top horizontal line) 
is obtained by the chain of dualities between various integrable systems. 
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bispectral duality, we shall connect the story to the Heisenberg SL{2) chain and to the 4d 
gauge theory]^] 

6.1 Liouville theory and rational Gaudin model 

Recall that the Liouville theory has central charge 

c=l + 6Q 2 , Q = b+l, (6.1) 

o 

and in the classical limit b — > oo so Q — > oo as well. Let us now consider a conformal 
block J 7 ^ a X ai (q) of the Virasoro algebra with central charge c — )■ oo with the four primary 
operators of dimensions 

Ai = a (Q-a ), A 2 = // (Q - A*o) , A 3 = ^(Q - fi{) , A 4 = ai(Q-ai), (6.2) 

inserted at points oo, 1, q, respectively on the S 2 with an intermediate s-channel state of 
dimension A = a(Q — a). In the above formula 

cto = lQ + Jlo, a = \Q + a, a 1 = lQ + Jl 1 , (6.3) 

where a is the SU(2) Coulomb branch coordinate. In the above formulae the mass parameters 
represent the following linear combinations of the SQCD quark masses mi^,3,4 

A*o = \{mi + m 2 ), Jio = |(mi - m 2 ), //i = §(m 3 + m 4 ), /Ii = §(m 3 - ra 4 ) . (6.4) 



There is an obvious notational conflict with [15], where /i's and m's are interchanged com- 
pared to our paper. We had to switch the notations in order to be consistent with Sec. |4j 
were m's are used for the quark masses. As far as the rest of the notations are concerned, 



they will be in agreement with 15 . Note that in Sec. |4j we treated all the four flavors as 
fundamental hypermultiplets, however, in [15] as well as in |11| two of them, with masses 
m 3 and m 4 are considered to be fundamental and two others, with masses m\ and m 2 to 
be antifundamental. For the purposes of Sec. [4] this turned out to be a mild difference and 
we were able to relate the 4d and 2d theories by studying the vortex effective theory. Also 
from the GLSM perspective it was natural to distinguish fundamental and antifundamental 
fields. In this section we have to be more careful about this issue as contributions from 
the fundamental and anti-fundamental multiplets to the Nekrasov partition at finite e are 
different. 



Note that all conformal dimensions (6.2) diverge at least linearly with b, however, as 
we shall later see, in order to match the Liouville CFT with the four dimensional theory 
in this limit, the dimensions will diverge quadratically and proper regularization is needed. 



Teschner in |62| have identified effective twisted superpotential (1.1) with the NS limit of a 



Liouville conformal block on the sphere as well as the proper regularization of the conformal 



9 Some remarks on the role of bispectral duality between the trigonometric Gaudin model and the XXX 
chain on the classical level can be found in [61] 

10 According to the NS dictionary this is also a Yang- Yang function 



35 



dimensions. Conformal block as the function of punctures' locations was found to 

satisfy the KZ equation (5.9) for the dual WZNW model with level k and b 2 = —(k + 2) _1 



1 d^izi) 

= T-LGaud^{zi) , i = l,...,L, (6.5) 



k + 2 d,Zi 



where "Hcaud is the Hamiltonian of the rational Gaudin modef^] Thus the large b limit 
corresponds to taking k — > —2. The conformal dimensions of chiral primary operators get 
rescaled and become 

* = (6-6) 



as b — > oo . For S 2 with four punctures at oo, 1, q and respectively from ( 6.2|6.3 ) and (6.17) 
we obtain 



x b 2J \ b 2 

as b — » oo. Our next step is to allow the mass parameters fi a and Jl a scale with b upon 
identification with the 4d theory. 

Equivalently one can also probe Liouville conformal blocks with surface operator inser- 
tions 



64 , those conformal blocks satisfy Gaudin eigenvalue problem in the NS limit 
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6.2 M = 2 SQCD in the NS Omega background 

On the 4d gauge theory side, we compute the Nekrasov partition function for the 4d M = 2 
SQCD with mass parameters //o, fio, Jj-i whose instanton part is 

iWa, Mo ,/Io,Mi,£i) = C 1 - ti^^-^rgPcM • (6- 8 ) 

where a = |Q — a and a is the SU{2) Coulomb modulus. For a generic Omega background 
the AGT dictionary says the deformation parameters are related to the 2d Liouville theory 
via 

b 2 = -, h 2 = ei e 2 . (6.9) 
£2 

The usual NS limit €2 — > corresponds to b — > 00 and e\ is kept fixed, then the Liouville 
theory becomes classical as h — > 0. In this section, we will be rather interested in the 
quantum regime of the Liouville theory so we shall allow e x — > 00 such that h is kept fixed. 



11 See also 63 where the matrix model approach to Hitchin systems in connection with the Liouville CFT 



was constructed. 



12 See also 65 where a systematic study of 4d gauge theories with surface operators in Omega-background 
was done. 
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It is clear that by a proper tuning of e\ and e 2 one can obtain any desired value of the Planck 
constant. 

As we have already discussed above, in the NS limit a more appropriate object to study 



is not the Nekrasov partition function but the effective twisted superpotential (1.1). As it 



was shown in [IT] that this superpotential also emerges from the (2, 2) GLSM which we have 
described in Sec. IH 

The DHL paper has done a perturbative calculation in the instanton number q in order to 



establish their 4d/2d duality (4.5) and the proof to all orders was further established in |28|. 



CDHL showed that in the NS limit the Nekrasov partition function can be represented as 
an integral over a finite set of variables and can be explicitly evaluated, and the saddle 
point condition is shown to be equivalent to the Bethe ansatz equations for the SL(2) XXX 
chain. One may ask immediately why the vortices are relevant, indeed they only exist in a 
Higgs branch of the four dimensional theory, whereas the AGT statement relates Liouville 
momenta with Coulomb branch coordinates. In order to understand this, let us recall that 
at zero value of the FI term the Higgs branch touches the Coulomb branch, and as it was 
by making a proper limit in the relatior^ ; 



pointed out in 



11 



m 2+a - n a e 



1,2. 



(6.10) 



one may recover any point of the Coulomb branch of the U(2) SQCD. Indeed, as e — > 
the Higgs lattice becomes more and more dense filling the Coulomb branch in that limit. 
However, for what we are doing here, the opposite e — > 00 limit is relevant, as it is required 
by the connection to the Liouville theory. Still we want to be able to cover any point on the 
Coulomb branch, so one has to scale the fundamental masses m a with e as well in order to 



keep combination (6.10) finite. So at any given Liouville momentum we only need to sit at 



a certain point on a Coulomb branch and the Higgs branch root has all information we need 
about that point. Recall that the anti-fundamental masses and, correspondingly fi and Jl 



are not affected by (6.10) and therefore do not scale with e. 



We now make an observation that the ground state equations for the (2,2) GLSM (5.1) 



or second equation in (5.8) where my are now denoted as my (antifundamental) and my 



became 777.3,4 (fundamental) respectively 



11 



n 

a=l 



\i - m 2+a + |e 
A,; - m n - \e 



K 

n 

3=1 



Xi — Ai 



Xj + e 



(6.1i; 



can be written as the second equation from the MTV dual pair (5.8). In order to see this we 



need to employ (6.10) and substitute 7773 and 7774 into the numerators of the let hand side of 



(6.11). Then we take the limit of large e keeping in mind that rapidities Xi also scale with e. 



Neither Coulomb moduli a a nor the antifundamental masses my enjoy this scaling, so they 
will drop out from the equations. We then arrive to (5.8) where z 2 jz\ = q and 



n 



a 



K a + 2. 



(6.12) 



13 From now on we shall work with the U(2) SQCD with 4 flavors. 

14 The obvious notational conflict occurs here. In spite of this we still keep the notations of this section 
and Sec. [4] as they are since both are natural in where they stand. We hope this issue will not confuse the 
reader. 



37 



6.3 The duality 



Now let us start connecting the story with the Liouville. By means of the bispectral duality 



these equations are mapped onto (5.4) yielding the trigonometric Gaudin model from the 



Heisenberg chain. Note that (5.4) depends only on two points Z\ and z 2 corresponding to 



the locations of the NS5 branes in 6-10 plane in Fig. |4| However, the Liouville conformal 
block depends on four operators sitting at oo, 1,5,0. Interestingly we can mention that 



trigonometric Gaudin Bethe equations (5.4) when only z\ and z 2 punctures are involved can 
be treated as a rational sl(2) Gaudin Bethe equations on S 2 with all four punctures included. 
Indeed, 

(6.13, 



6=0 



Zb 



where % 



1,2,3 



{oo, 1, q, 0}, is equivalent to (5.4) with 



eu 2 = K, ez/ 3 = m 3 — m 4 — e = 2/^ — e 



(6.14) 



being spins of the s[(2) representations sitting at points q and 0. Specification of z/ is not 
important, as the corresponding contribution drops out from the equation since zq = oo. 



Also, as we have already mentioned in (6.12) there is an exact matching between the number 
of the Gaudin Bethe roots with the parameters n a of the Higgs branch. In other words, all 
sectors of the trigonometric Gaudin model Hilbert space parameterized by number of Bethe 
roots (excitations over the Bethe vacuum), by means of the bispectral duality, are mapped 
onto various points of the Higgs branch lattice {n a } of the four dimensional theory. Finally, 



the value of v\ can be found from (5.6) and (6.12). We conclude that v\ = —2, which formally 
corresponds to the spin —1 representation for z\ = \ 

Note that one should also take out the U(l) factor from the U{2) gauge group, as it does 
not have an analogue in the Liouville theory. Imposing it on the U(2) Coulomb moduli a 1; a 2 



with the help of (6.10) we get 



m 3 + m 4 - (m + n 2 )e = 



or, using the Liouville mass parameters, one gets 



Hi _ n\ + n 2 
T ~ 2 



(6.15) 



(6.16) 



The U(l) condition balances the count of the parameters on both sides of the correspondence 
as in order to match st(2) spin at z 4 = we used only one antifundamental mass parameter 
(which is related to the fundamental one). 

Again, from the gauge theory perspective we are interested in keeping Coulomb branch 



parameters in (6.10) finite while masses /iq and fii and e are sent to infinity. The rescaled 



conformal dimensions (6.7) upon identification b = e and by using (6.16) therefore read 



5n = — 



8i=0, 5 2 



K (K 



+ 1 



#3 = 7o(7o + 1) 



(6.17) 
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7o 



(6.18) 



where 

nx fiv _ 1 _ U3 
" 2 2 2 ~ 2 ' 

and we recall K — h\ + h 2 is the number of the D2 branes stretched between the NS5 
brane at z 2 = q and the D4 branes. We can also see that in such limit the SU(2) Coulomb 
coordinate and the anti-fundamental masses have dropped from the formulae. In the last 



two terms of (6.17) we recognize sl(2) quadratic Casimir eigenvalues on representations of 



spins and y respectively. We can see from (6.14) and ( |6.17 ) that the matching occurs 



at these points. Vanishing eigenvalues 5± confirms the fact that v\ = —2 corresponds to the 
spin —1 representation. Spin of the representation at z = 00 is formally equal to — |. 

Here is the summary table of the correspondence between the objects we have discussed 
in this section in addition to the standard AGT dictionary 



Liouville conformal block at b — > 00 
on S 2 with four punctures 


U{2) , N f = 4 SQCD instanton 
partition function in the NS limit 


Rational Gaudin model from KZ 
equation on conformal blocks 


SL(2) spin chain from the ground state 
equation for the 2d GLSM dual to 4d theory 


Punctures' positions z 2 /zi 


Instanton number q 


sl(2) spin at z 2 = q 


U(l) condition, number of D2 branes 
emerging at NS5 brane at z 2 = q 


Conformal dimensions of chiral operators 
at points z 2 = q, z 3 = 


Quadratic sl(2) Casimir eigenvalues on 
spin ^hi + \fi 2 and 
— \h\ + \fi 2 — | representations 


Gaudin Hilbert space sectors with 
different number K a of Bethe roots 


Higgs branch lattice {n a } 



6.4 Generalization to SU(2) linear quivers 

One can easily generalize the above construction to the Liouville theory on S* 2 with L + 3 
punctures. A natural quiver gauge theory associated to this Riemann surface has L SU (2) 
gauge nodes with Coulomb moduli ai successively connected together. Liouville conformal 
block of L + 3 operators located at points 

00, 1, q u qiq 2 , qiq 2 ■ ■ ■ qL, , (6.19) 

with the following scaling dimensions 

a (Q-a ), fj. (Q - fa) , . . . , m L (Q-m L ), a L+1 (Q - a L+1 ) , (6.20) 

respectively glued by operators of dimensions a>i(Q — «j) in the intermediate s-channels. 

NS limit of such quiver theory has been elaborated in [28]. Brane interpretation of the 
bispectral duality is very useful in this case. The corresponding Hanany-Witten picture view 
from "below" (in (6 + zlO) — 7 space) is shown in Fig. 10 Quiver theories have bifundamental 



matter with masses //^ , so the Higgs branch conditions get changed 



,(p) 



(P) 



(6.21) 



k=l 
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2D4'.s 1 2D4's 2D4's 



2D4',s 



2 D4's 2 D4's 



OG 



NS5 



NS5i 



A'i _D2's 



21 



NS5 



if 2 £>2's 
22 NS5 



A" 3 D2's 

^3 





K L D2's 
NS5 z L 



Figure 10: (6 + HQ) — 7 slice of the CDHL quiver construction. Vertical lines correspond to stacks 
of D2 branes. 



Using the above relation we can express conformal dimensions (6.20) in terms of Coulomb 
branch coordinates, quantization parameters and bifundamental masses. Performing the 



rescaling analogous to (6.7) we conclude that operators located at 



zi = qi, z 2 = qiq 2 , z 3 = qiq 2 q^ 



have dimensions 



+ 1 



Kl 
2 



K T 



z L = qi...q L 



+ 1 



(6.22) 



(6.23) 



where K,, 



(0 , aW 



n\ +n 



- 1, . . . , L, corresponding to the sl(2) Casimir eigenvalues on repre- 
sentations of spins \K\ . . . , \Kl- Spins sitting at each point Z{ (multiplied by 2) correspond 
to the total number of D2 branes stretched between the z'th NS5 brane and both D4 branes. 



Our construction has to be supplemented by L £7(1) conditions similar to (6.16) for each 
gauge group. 

The full treatment of the AGT duality for linear quivers requires the construction of the 
bispectral dual to the twisted anisotropic SL(L,M) chain, which emerges from it 28 . We 
shall postpone this analysis for the future work. 



7 Conclusions and Outlook 

In this paper we have investigated BPS solitons - strings, monopoles and domain walls in 
M = 2 four dimensional gauge theories in Omega background with the Nekrasov-Shatashvili 
limit imposed. We derived the central charges for these solitons from the supersymmetry 
algebra and observed that string and domain wall charges are proportional to the external 
graviphoton field present in Omega background. At large values of e, string and domain wall 
tensions are large which makes semiclassical considerations legitimate. Existence of a BPS 
string in pure SYM implies fractional windings for the adjoint scalar and hence presence of 
a conical singularity. We have presented arguments that its tension is nevertheless finite. 



Moreover we have argued that BPS monopoles found in 14 are actually located on domain 



walls interpolating between two different vacua of the theory. The next step along this road 
will be to analyze the complete set of the moduli space of the corresponding solutions and 
worldvolume theories of BPS solitons. Also one may expect more surprises in study of wall 
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crossing phenomena for such BPS objects. We postpone the discussion on these issues for a 
separate publication. 

There is an interesting question regarding the behavior of solitonic BPS states at small 
e. In the pure SYM strings and domain walls become highly quantum objects and could 
potentially condense. However it is not clear if such condensation of the extended defects 
should be taken into account. In SQCD at large value of the FI term in the large e limit the 
standard nonabelian string is recovered. Note that small e limit corresponds to semiclassical 
regime from the integrability viewpoint, hence the potential condensation of BPS solitons 
also deserves further investigation in the Hamilton- Jacobi framework. 

An interplay between quantum integrability and Omega deformed gauged theories allowed 
us to apply some known dualities in the integrable systems literature. Degrees of freedom 
in an integrable systems are interpreted as coordinates of the corresponding brane positions 
since all dualities can be reformulated in brane language. The integrability being just the 
reflection of the symmetry of the brane geometry involved plays a role of a consistency 
condition for the whole construction. We have shown that different dualities connecting 
Gaudin, XXX and Calogero type models can be explained in terms of brane geometry both 
at classical and at quantum levels. 

As a byproduct of the vortex construction and the dualities between integrable systems 
we were able to reconstruct the AGT correspondence in the NS limit. The Liouville CFT 
has infinite central charge and all scaling dimensions needed to be regularized. In the paper 
we discussed the Liouville theory on S 2 with L + 3 punctures, which was shown to be dual to 
linear quiver gauge theories with corresponding matter content. Our construction certainly 
has to be extended to other known AGT dual pairs, e.g. a torus with multiple punctures, 
etc. 

There are many questions which certainly deserve additional study. Surely, more compli- 
cated defects involving strings, domain walls and monopoles have to be explored. It would 
be interesting to investigate similar defects in five and six dimensions and with the complete 
Omega deformation beyond the NS limit. 
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A The Gaudin Model 



Here we discuss the Gaudin model - the key tool in our AGT construction, its relations with 
the XXX spin chain. 

Gaudin model from XXX chain. The Gaudin model is the simplest example of the 



Hitchin system on a sphere with marked points 66 . It is also known to be a large impurity 
limit of an anisotropic twisted XXX spin chain. This fact can be realized both in the transfer 
matrix at the classical limit and in the Bethe ansatz equations in the quantum case. We 
shall be interested in the quantum case and upon the proper limit Bethe ansatz equations 
for the Gaudin model can be obtained. Let us start with Bethe equations for anisotropic 



XXX s spin chair 15 with twist q = e 



a=l z j=l J 

By taking logarithms of both parts of the above equations, then rescaling 

t 

Xi i y xXi, v a i y xu a , f - , (A. 2) 

x 

and sending x — > oo we arrive at the following set of equations 

1 N c K o 

a=l 3=1 J 

which are nothing but Bethe equations for the Gaudin model. The anisotropies u a at each 
site still play the role of the inhomogenities in the model, while the twist q in the XXX chain 
play the role of the external field in the Gaudin system. As we can see the latter vanishes 
as e — > oo. 

Bethe ansatz equations for the Gaudin model. Let us now recall how the Bethe 
ansatz equations for the rational Gaudin model with the Lie algebra symmetry q are derived. 
For our purposes we need merely $j = s((2) and and L points on the sphere. At each 
point we fix a representation V^(Vi), . . . , V{v£) of s((2) algebra with some dominant weights 
v a , a = 1, . . . , L. According to the Bethe ansatz prescription |36| we construct the following 
operator 

0=1 " a a=l \ U Zal 

where l-i a are Gaudin Hamiltonians at each site of the lattice 

dim(fl) (6 ) (6) 

«. = EEfV- (A - 5) 

b=£a a=l 



lo 



We measure spectral parameters in units of it here. 
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where 04 of the acts with $ a G s((2) on the 6-th site of the spin chain and with identity 
on the others. A(V a ) are eigenvalue of the U(s((2)) quadratic Casimir acting on V(v a ). For 
such a system Bethe ansatz equations for the sector with n a Bethe roots read as follows 

6=1 1 b j=l 1 3 

B Supersymmetry Algebra and Central Charges 

M = 2 supersymmetry algebra in four dimensions has the following form 



} J p} — £a/3t IJ Z mon + (Zd.w. ',, , 



{QlQi} = + (z d (b.i; 



There are three types on central charges: string, monopole and domain wall types. The full 
global symmetry of the theory is SU(2) L x SU{2) R x SU(2) n x SU(2) C . It is broken by the 
Omega background in the NS limit to SU{2)r + h x SU(2) c . Twisted supercharges 

Q = 5j Qq, , Qm = ipm) Qla ; Qmn = i&mn) iQ a • (^-2) 

The former operator above is also known as BRST operator. The transformations can be 
inverted as 



Plugging these formulae into (B.I ) we get the twisted version of the supersymmetry algebra. 



C Some Notations 

In this work we are intended to use Euclidean signature. We benefit from this while studying 
static configurations, where in the gauge A4 = the Lagrangian is nothing but the energy 
density. Below we list some definitions and conventions. Raising and lowering of spinor 
indices is performed by means of Levi-Civita symbol 

e 12 = -e 21 = -e 12 = e 2 i = l, (C.l) 

the same for the Levi-Civita symbol with dotted indices. The definition is consistent with 

e a/3 e/3 7 = ^- (C.2) 

Scalars can then be obtained by contracting spinor indices, for example r) a x a = ~V a Xa- 
Vector indices are contracted with Euclidean metric 5 mn . 
Sigma matrices 

(<T m )aa = {-tr\ -IT\ ~%t\ 1) , 

(v m h* = (ir\ it\ it\ 1) = ((a m ) aA y , (C.3) 
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where r 1,2 ' 3 are standard Pauli matrices. Thus for a m the undotted index goes first, whereas 
for a m it is the last. 
Lorentz projectors 



/ mn 



1 / 1 _m\ l — n\a l _nN /-m\ci \ 

(* mn U = \ ({v m )a° n Yp - (v n )a° m r $ ) • (c.4) 

Chiral and antichiral electromagnetic field strength 

F afi = -;F„(ff m ) m (/); = -{E + B) ■ a a 2 , 

= \F mn {a m f a {a n f a = (-E + B) ■ a a 2 . (C.5) 
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